Gurevich-Zybin system 



Maxim V.Pavlov 

Abstract. We present three different linearizable extensions of tfie Gurevicli-Zybin 
system. Their general solutions are found by reciprocal transformations. 

In this paper we rewrite the Gurevich-Zybin system as a Monge- Ampere equation. 
By application of reciprocal transformation this equation is linearized. Infinitely many 
local Hamiltonian structures, local Lagrangian representations, local conservation laws 
and local commuting flows are found. Moreover, all commuting flows can be written as 
Monge- Ampere equations similar to the Gurevich-Zybin system. 

The Gurevich-Zybin system describes the formation of a large scale structures in 
the Universe. The second harmonic wave generation is known in nonlinear optics. In 
this paper we prove that the Gurevich-Zybin system is equivalent to a degenerate case 
of the second harmonic generation. Thus, the Gurevich-Zybin system is recognized as a 
degenerate first negative flow of two-component Harry Dym hierarchy up to two Miura 
type transformations. A reciprocal transformation between the Gurevich-Zybin system 
and degenerate case of the second harmonic generation system is found. A new solution 
for the second harmonic generation is presented in implicit form. 
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1. Introduction 



The invisible nondissipative dark matter plays a decisive role in the formation of a 
large scale structures in the Universe: galaxies, clusters of galaxies, super clusters. Corre- 
sponding nonlinear dynamics can be described (see [14]) by following hydro dynamic-like 
system 

p^ + V(pu) = 0, ut + (uV)m + V$ = 0, A<l> = p, (1) 

where the first two equations are usual hydrodynamic equations (the continuity equation 
and the Euler equation, respectively), but the third equation is the famous Poisson equa- 
tion. At least for the first time this system was derived by J.H. Jeans (see [16] and also 
[31]) for a description of instabilities of a homogeneous distribution of a matter. 

Such dynamics of dissipationless gravitating gas is a special limit {e 0) of another 
system (another sign of p is unessential) 

p^ + V(pu) = 0, ut + (uV)m + V$ = 0, A$ = ee* - p 

describing fully nonlinear flows in a two temperature unmagnetized coUisionless plasma 
in dimensionless variables (nonlinear ion-acoustic waves; see for instance [30]). 

The main advantage of the Jeans theory is a reckoning of two factors: a gravity 
attracting a matter in separate lumps and clots; a pressure decreasing an inhomogeneity 
of a matter in the Universe. 

Recently a new achievement in the investigation of the system (1) was made (see [14]) 
in Cosmology. The nonlinear one-dimensional dynamics of a dark matter is described by 
the equations [14] 

Ut -\- uUx -|- 1; = 0, Vt + uVx = 0, (2) 

where p = Vx-, v = ^x- The analysis of equations (2) in a multimode form demonstrates 
the transition from the hydrodynamic to the equilibrium kinetic state [14]. It means, 
that the exact solution of the equations (2) describes a fundamental physical process (see 
[14] again). 

It is amazing that the inhomogeneous hydrodynamic type system (2) can be 
integrated, up to the first singularity, by the Hodograph Method (see [14]). For this 
reason further we will call the system (2) the Gurevich-Zybin system emphasizing that 
the one-dimensional reduction (2) of the system (1) is integrable. 

Here we give the general solution by the method of Reciprocal Transformations. More- 
over, we present three different linearizable extensions of this system (2) with their gen- 
eral solutions given by corresponding reciprocal transformations. Actually these recip- 
rocal transformations have clear pure mathematical (hodograph method) and physical 
(transition from Euler to Lagrange variables) interpretations. In the next section we 
present three linearizable extensions of the Gurevich-Zybin systems with their general 
solutions. In second section relationship between two-component generalization of the 
Hunter-Saxton equation and the Gurevich-Zybin system is established. In the third sec- 
tion the Gurevich-Zybin system is rewritten as a Monge- Ampere equation (following the 
approach developed by Andrea Donato). In the fourth section a bi-Hamiltonian struc- 
ture of the Gurevich-Zybin system is found (following the approach developed by Yavuz 
Nutku). In the fifth section, by the application of a reciprocal transformation, a simplest 
recursion operator is constructed. Infinitely many local conservation laws, local commut- 
ing flows, local Lagrangians and local Hamiltonians are found. Moreover, all commuting 
flows are Monge- Ampere equations. Thus, the Gurevich-Zybin system is a member of an 



integrable hierarchy of Monge- Ampere equations. In sixth section a bi-Hamiltonian for- 
mulation for the Gurevich-Zybin system is given in a canonical form. The Gurevich-Zybin 
system is recognized as a first negative flow of two-component Harry Dym hierarchy. In 
seventh section Miura type and reciprocal transformations between the Gurevich-Zybin 
system and Kaup-Boussinesq hierarchy are given. In eighth section we finally prove that 
the Gurevich-Zybin system is equivalent to a degenerate case of the second harmonic 
generation system up to above-mentioned transformations. A new solutions of the second 
harmonic generation system is found. In conclusion we discuss about sort of integrable 
problems belonging some different hierarchies of integrable equations. 

2. General Solution 

The Gurevich-Zybin system (1) in one dimensional case precisely has a form 

Pt + dx{pu) = 0, ut + uU;, + ^^ = 0, ^xx = P- (3) 

This system can be generalized at least in three different forms: 

Pt + d^{pu) = 0, ut + uu^ + p'"{(!>^) = 0, $xx=P, (4) 
Pi + d^{pu) = 0, ut + uu^ + <l>^ = 0, <^^^ = B{p), (5) 

p, + d,{pu) = 0, + + = 0, d,C{<^,)=p, (6) 

where p{z), B{p) and C{z) are arbitrary functions. It is interesting that is not obvious 
that the system (4) is equivalent to the Gurevich-Zybin system (3). Indeed, the system 
(4) written hke (2) 

Zt + uZx = 0, Ut + uUx + p"'{z) = (7) 

is exactly (2) up to point transformation v = p"\z). One can introduce reciprocal 
transformation 

dz = pdx — pudt, dr = dt. (8) 
Then dx = pdz and dt = dr — pudz- Thus, the system (4) has the form 

(^)r = Ur = -p"'{z), (9) 

in new variables, where z = Thus, the general solution of the system (9) is 

u = -p"\z)r + D\z), - = -p""{zy/2 + D"{z)r + E"iz), 

P 

where D{z) and E{z) are arbitrary functions. Finally the general solution of the system 
(4) can be given in implicit form 

u = -p"\z)t + D'{z), p=[-p""{z)ty2 + D"{z)t + E"{z)]-\ 

X = -p"'{z)e/2 + D'{z)t + E'{z), (10) 

$ = {p"{z)-zp"'{z))t'l2 + {zD\z)-D{z))t + zE\z)-E{z), 

where 2; is a parameter here. 

Above reciprocal transformation applied to the system (5) yields 

(-)t = Ur = -V, Vz = B{p)/p, 

P 



where v = If function p can be explicitly expressed from algebraic equation 

p 

d9 



G{z) 



9\ 2j B{9)9-^d9- F{ 



where F{z) and G{z) are arbitrary functions, then the general solution of the system (5) 
can be obtained. For instance, if B{p) = p, then a general solution is already given by 
(10) (remember that in such sub-case systems (4) and (5) coincide if p{z) = z^/2A); in 
simplest perturbed case B{p) = p + 6/p {6 = const) the general solution is expressed via 
Weiershtrass elliptic functions 

P = -G{z), —F{z))\ , M^ = -p, v^ = \-—<p. 
The reciprocal transformation (8) applied to the system (6) yields 

(-)r = Mz, Ur = -p"'{v), C{v) = Z, 

P 

where v = Thus, the general solution of this system is 

V = V{z), u = D'{z) - I2"'{v)t, 



p = [E"{z)+D"{z)t- 



p"'\v) 



-2l-l 



r 



2C'{v 

X = E'{z)+D'{z)T-i2"'{vy/2, 

where D{z), E{z) are arbitrary functions and V{z) is inverse function to C{v). Finally 
the general solution of the system (6) can be given in implicit form 

V = V{z), u = D'{z) - p"'{v)t, 

p = [E'\z)+D"iz)t-Q^^tr\ 
X = E\z)+D\z)t- p"\v)t'l2, 

$ = zE'{z)-E{z) + {zD'{z)-D{z))t+[G{v)-zG'{v)]e/2, 

where G' = p"'C' and z is a parameter here. Also the system (6) can be written in 
hydro dynamic like form 

p^ + ^x{pu)=Q, Ut + uUx + ^dxP = 0, (11) 

where the pressure P is a nonlocal function of the density p 

p = p{v{d~'p)). 

In particular case C{v) = v the system (6) coincides with the system (4); the system (6) 
was written in form (11) in [6] for particular case p{v) = t;^/24 and C{v) = v. 



3. Two-component generalization of the Calogero equation 

The simplest two-component linearizable generalization 

Vt + dxivu) = 0, u^t + uu^x + ^(?7, Ux) = (12) 
of the Calogero equation (see [7]) 

Uxt = UUxx + R{Ux) 

was presented (functions '^{a^h) and R{c) are arbitrary here) in [26]. Some particular 
cases of the Calogero equation like the Hunter-Saxton equation (i.e. i?(c) = vc^, where 
v = const; see [15]) are interested from physical point of view. The Calogero equation 
is linearizable by a reciprocal transformation (see [26]). For instance, the Hunter-Saxton 
equation is related to the famous Liouville equation by a reciprocal transformation (see 
[8]). Thus, the system (12) is a natural generalization of the Liouville equation on two- 
component case up to module of a (invertible) reciprocal transformation. 
The reciprocal transformation 

dC, = rjdx — r]udt, dy = dt (13) 
applied to the system (12) yields the ordinary differential equation 

^yy ~^ ^(c , Sy) 0, 
where s = — ln?7, which can be reduced to the first order equation 

arfa + ^'(e"^ a)rfs = 0, (14) 
where a = Sy = Ux- Then a general solution can be constructed in two steps from 

uc = (-)y, dx = -d( + udy. 
rj rj 



(15) 



A solution q{s, a) of the linear equation 

dq \E'(e^'*, a) dq 
ds a da 

determined by the characteristic equation (14) yields extra conservation law 

Pt + dx{pu) =0, 

where p(?7, Ux) = rjexpq. The comparison of the second equation in (12) and the second 
equation in (3) yields another relationship 

Thus, a solution of the nonlinear equation (substitute \1/ from above equation into (15)) 

qs = {a-\ )qa 

a 

describes a transformation between the Gurevich-Zybin system (3) and two-component 
generalization of the Hunter-Saxton equation (12). 
Remark 1: Above equation under the substitution 

q = s + ln(^e"^' - a^) 



transforms into well-known inhomogenious Riemann-Monge-Hopf equation 

2c 



Uy + nUc 



v 

where y = e~^*/3 and c = . Its general solution can be given just in parametric form 

where ^ is a parameter and Ai{C,), ^2(0 ^'"^ arbitrary functions. However, the general 
solution of the equation of the first order depends by one function of a single variable 
only. It means, that if for instance Ai{^) ^ const, then by re-scaling ^ the 

general solution take a form 

n = l^y"'/' + lAiOy-'/^ c = iy'l' + A{i)y''\ 

where A{^) is an arbitrary function. 

In a particular case the substitutions (see [26]) 

P = -\{nl + r/2), ^{71, u,) = \{ul - r/2) (16) 

connect the Gurevich-Zybin system (3) with two-component Hunter-Saxton system (see 
(12) and [25]), which is a bi-directional version of the Hunter-Saxton equation. A general 
solution of the Gurevich-Zybin system in field variables rj and u has the implicit form 
with respect to the parameter C, 

^ = [^+i^'(o(t-m(c))rs 

u = ItkiO -\j miOdkiO, (17) 

X = /[^ + \k'iOm\0]dC + \t'k{C) -\tj m{C)dk{C), 

where m{(), k{Q are arbitrary functions. One can substitute above expression for t] to 
the first equation in (16) and compare expressions for (p, x, u) from (10) (in this case 
ti"'{z) = z) and above, then a relationship between the arbitary functions /c, m and E 
will be reconstructed. 

4. Reformulation of the Gurevich-Zybin system as a Monge- Ampere 

equation 

A lot of physically motivated nonlinear systems can be written as Monge-Ampere 
equations (see [29]). To this moment we have no unique method for constructing such 
relationships. One simple approach was suggested to author by Andrea Donato (see [9]) 
at the "Lie Group Analysis" conference in Johannesburg (South Africa) at 1996. 

The Gurevich-Zybin system in physical field variables (3) has four local conservation 
laws 



Ut + d^{u'/2 + <^) = 0, p, + d,{pu) = 0, 
dt{pu) + d^{pu' + ^l/2) = 0, dtipu^ - ^l) + d^ipu') = 



as a consequence of obvious local Hamiltonian structure 

where the Hamiltonian is H2 = ^ Jl—u^Vx + v'^]dx, the momentum is ifi = / uvxdx, 
and two Casimirs are functionals Qi = J udx and Q2 = J pdx of corresponding Poisson 
bracket 

{p(x), u{x')} = {u{x), p{x')} = 5'{x - x'). (19) 

The existence of above first three local conservation laws is obvious. However, the fourth 
conservation law is not easy to check. Since p = ^xx, then pu = —^xt, then pu^ + $^,/2 = 
then above fourth conservation law is valid. Eliminating physical field variables p 
and u from these three equations, the Monge-Ampere equation is given by 

- = Ul<^>xx. (20) 

In paper [14] the Gurevich-Zybin system was linearized by a hodograph method. A 
general solution has been presented too. Thus, this Monge-Ampere equation is linearizable 
and has the general solution in implicit form (see the end of section 2) 

$ = -^vH^ + {vD'{v) - D{v))t + vE'{v) - E{v), x = -vt^/2 + D'{v)t + E\v), 

where D{v) and E{v) are arbitrary functions, t; is a parameter here. 

Since the Gurevich-Zybin system can be written in the form (4), (7), we shall use an 
arbitrary value of function p{z) in next two sections. 

5. Bi-Hamiltonian structure 

The Gurevich-Zybin system (7) has local bi-Hamiltonian structure, where the first 
local Hamiltonian structure is (18) 

m m .... 

just three conservation laws are 

dtp + dxipu) = 0, 

dt{pu) + dx[pu^ + p"{^x)] = 0, 

dt[pu^-2fi"{^x)] + dxW) = 0, 

associated with the first Poisson bracket (19). The Hamiltonian is H2 = /[— pw^ + 
2n"{^x)]dx, the momentum is i^i = J pudx and the Casimir is Hq = J pdx. Correspond- 
ing Lagrangian representation is 

•^1 = P"{z)]dxdt, (22) 

where u = —Zt/zx- Thus, the Lagrangian 

creates the Euler-Lagrange equation 

<!>xx^tt - ^It = /i"($x)$xx, (23) 



which is a Monge- Ampere equation (cf. (20)). 

At the same time (23) allows other Lagrangian representation (see [24]) 

^2 = / [l^..^t - l^{^x)]dxdt. (24) 



2 

Thus, the Monge-Ampere equation (23) has the second Hamiltonian structure 

ft — dx—^ ; Zt — Ox— 

oz or 
determined by the local Poisson bracket 

{z{x), r{x')}2 = {r(x), z{x')}2 = 5'{x - x'), (25) 

where r = ^^x^t, the Hamiltonian is H2 = J[-^+fi{z)]dx, the momentum is Hi = J rzdx, 
two Casimirs are Qi = J rdx and Q2 = J zdx. Four local conservation laws associate with 
above second Hamiltonian structure are 

2 -1 3 

zt = dx{-), d,[^ + ^i{z)]=dxW{z)- + -dx{^)], 

Zx Zx Zr^ 

n = dx[fx'{z) + dxi-^)], dt{rz) = dx[zfx'{z) - fx{z) + -dx{—)]. 

^Zrr< ^ Z, 



X 



6. Recursion Operator. Integrability of the GZ hierarchy 

Applying the reciprocal transformation (8) simultaneously to both Lagrangian rep- 
resentations (22) and (24), one obtains the variation principles in another independent 
variables 

Si = j [\xl + ^^"{z)x,]dzdT and 82 = j {U>1 - ^(z)$..]rfzrfr, 

where = (see the first conservation law associated with the second 

Hamiltonian structure (25): d^ = zdx + ^dt or c/$ = xdz — ^dr, where xz = $ + <!>). Then 
the Euler-Lagrange equation is Xrr = ~l^"'{.z). This equation can be easily integrated (see 
(10)). Corresponding Poisson brackets (see (19) and (25)) 

{x{z), u{z')}i = -{u{z), x{z')}i = 5{z - z'), 

{piz),^z')}2 = -{^z),p{z')}2 = 5iz-z'), 

where u = and (i$ = xdz — pdr (i.e. p = r/p = $(), create the recursion operator 

5? 



R 
where 



{x{z), u{z')}2 = -Mz), x{z')}2 = -5"{z - z'). 

Thus, the Gurevich-Zybin system in these independent variables has an infinite set of 
local Hamiltonian structures, conservation laws and commuting flows. For instance, all 
such Hamiltonians are 

Hk = {-If j [^p^')' + p^''^^\z)¥'''^]dz, A; = 0, ±1, ±2, ... 



Corresponding commuting flows 

can be easily integrated (see (10)). However, in the independent variables {x, t^) they 
can be written in the form (cf. (7)) 

dtku'' + u^d^u^ + /x(2'=+^)(z) = 0, dt^z + u^d^z = 0, (26) 

where 

mO = m, u^+^ = -d^u'', u-^-^ = d-^{pu-^), A; = 0, 1, 2, ... (27) 
P 

Thus, all commuting flows to the Gurevich-Zybin system created by above bi-Hamiltonian 
structure are Monge-Ampere equations (cf. (23)) 

where 

= $^tfe = k = 0, ±1, ±2, ... 

All local Lagrangians are 



S2,k = J[^^,.<^l-fl^''\<^.)]dxdt' 

Si,k = J[^+fi^"'^'\<^.)]dxdt^ 

.' ^^xx ^xx ^xx 



Remark 2: All commuting flows have infinitely many different local representations via 
different pairs of field variables {z, u^), see (27). For instance (cf. (26)) 

zt. + d,u^"' = 0, dt.u''-' + ^^^^ + fi^^'^'\z) = 0, 

^x 

zt^ + d^i^) = 0, d,.u'-' + + fi^''^'\z) = 0. 

Zx IZ,^ 

The theory of integrable systems with a multi-Lagrangian structure is presented in [23]. 
Usually, every local Lagrangian creates a nonlocal Hamiltonian structure. Such explicit 
formulas of nonlocal Hamiltonian structures, nonlocal commuting flows, nonlocal conser- 
vation laws as well as nonlocal Lagrangians can be found iteratively from already given 
above formulas. 



7. Another bi-Hamiltonian structure 

Now in this and two next sections we identify f = z, i.e. we concentrate attention on 
case ii"'{z) = z, see (4)). In two previous sections we discuss bi-Hamiltonian structure of 
the Gurevich-Zybin system. Here we preserve the first Hamiltonian structure (see (18), 
(21)), but change the second one! In above section we proved that the Gurevich-Zybin 
system has infinitely many local Hamiltonian structures and Lagrangian representations (a 



general theory is presented in [23], see also [22]). However, this new second Hamiltonian 
structure (see below) is not from this set! 

The Gurevich-Zybin system (2) is an Euler-Lagrange equation of corresponding vari- 
ational principle (see (22), when ij!'{z) = /2) 

v2 



S = - f[^ + z^]dxdt. 
2 J Zx 



However, the astonished fact is that the Gurevich-Zybin system (2) has another 
Hamiltonian structure connected with the same Lagrangian density. Namely (see for 
details [22], especially formulas (43), (52-54) therein), the Lagrangian (cf. with S) 



S = l ! p^(2zt - px) + z']dxdt 

^ J Zx 



determines the same Euler-Lagrange equations (2) but with another Hamiltonian struc- 
ture 

where u = —px/zx (i.e. p = $t). 

Remark 3: This bi-Hamiltonian structure at first was discovered by Yavuz Nutku 
[21] and later it was independently found in [6] (see formula (9) therein) exactly as it 
was done in [22]. However, here we repeat and emphasize the main observation of this 
section is that both Hamiltonian structures have the same Lagrangian density] This is 
the first such example in the theory of integrable systems. 

Canonical representation for both Hamiltonian structures and recursion 
operator. The Poisson bracket 

{u{x),u{x)]i = 0, {p{x), u{x')}i = 6'{x - x'), 
{n(x), p(x')}i = 5'{x-x'), {p(x), p(x')}i=0 

of the first Hamiltonian structure is given in its canonical form (more details see in the 
review [10]). However, the Poisson bracket 

{u{x), u{x')}2 = —d^^6{x — x'), {u{x), p{x')}2 = —UxS{x — x'), 
{p{x), u{x')}2 = UxS{x - x'), {p(x), p{x')}2 = -{pdx + dxp)5{x - x') 

of the second Hamiltonian structure can be reduced by the Darboux theorem to the 
canonical form 

{w{x), w{x')}2 = 5'{x-x'), {ri{x), w{x')}2 = 0, 
{w{x), r]{x')}2 = 0, ivi^), Vix')}2 = S'{x - x') 
by the Miura type transformation (see the first equation in (16)) 

w = Ux, p = -^{w^ + r]^) . 

Then the Gurevich-Zybin system written in a modified form (see formula (26) in [26], 
other details in the last section 5 and references [6], [7] therein) 

r]^ + dx{ur]) = 0, Uxt + uUxx + = ^r]^ (28) 



can be recognized as the two-component generalization of the Hunter-Saxton equation 
(cf. [15], [8], [25]). 

Remark 4: In fact, the Casimir density t] of the second Hamihonian structure was 
found in [G] (see formula (18) therein). However, the Gurevich-Zybin system was not 
presented in form (28) there. Moreover, we emphasize the main result of this paper 
is that the Gurevich-Zybin system belongs to the well-known class of integrable systems. 
In this section we prove that the Gurevich-Zybin system is a member of an integrable 
hierarchy embedded into 2x2 spectral transform. 

Since the first Poisson bracket in new field variables has a form 

1 1 

{w{x),ri{x')}i = —[—6{x — x')]", {ri{x),w{x')}i = —6"{x~x' 



XII I 



{w(x),w(x')}i = 0, {77(x),?7(x')}i = -[^9^ + 9^^]5(x - x'), 

then the modified Gurevich-Zybin system (28) as a member of integrable hierarchy with 
all other commuting flows together can be written in bi-Hamiltonian form 

Vf^ = d^—. — = -dx-. 2—9^ + — xH— , 

07] 7] OW 7]'^ Tj'^ Of) 

An eigenvalue problem for the recursion operator as a ratio of both Hamiltonian structures 



i 



) = 2XdA ^1 



can be written as one equation 

if,,, + 4(A'r/2 + Aa)(/., + 2{\^rj^ + \a),^ = 0, 

where (p^ = ip,, = —2XrjLp and a = w,. However, above equation can be reduced to 

^xx + + Aa)V^ = 0, 

where (p = is a squared eigenfunction and ip, ip^ are linear conjugate solutions with 
different asymptotics at infinity A — * 00. This linear spectral problem (more precisely, 
just "a:"-dynamics) is well known in the theory of integrable systems: corresponding 
systems are members (commuting flows) of the two-component Harry Dym hierarchy (see 
for instance [2]). All such members of this hierarchy can be determined by the spectral 
transform 

^xx = -(AV + Afx)^, ipt = b^x - i^x^. (29) 

where 6(C, 77, A) is a 'polynomial function with respect to the spectral parameter A for 
positive members. The compatibility condition {4',,)t = ijPtlxx yields relationship 

(A^r/^ + \a), = l^dl + 2(AV + Aa)9,. + {XW + Aa) J6, 
where the two-component Harry Dym system (see [2]) 

Vti {~^)xy (^ti {~)xxx 



can be obtained if 6 = 2A/?7. Thus, twice potential two-component Harry Dym system 

Vt^ = ut, = (-)x (30) 

is the first member of positive part of above hierarchy and the first member of its negative 
part is the modified Gurevich-Zybin system (28) 

Vt-i + dx{ur]) = 0, u^t-i + uu^x + ^ul = \f (31) 

determined by the choice h = (2A)~^ — u (it means that we must identify t = t_i for the 
Gurevich-Zybin system (2)). 

Remark 5: The reciprocal transformation (see (8)) 

u 

dri = dti, dT_i = dt-i, dz = pdx — pudt_i — (—)xdti 

V 

simuhaneously hnearizes the Gurevich-Zybin system (see (9) and (10)) and preserves 
two-component Harry Dym system: 

w 1 _ w _ 1 _ 

Pti = -{-)xx, Wt, = {-)xx, Vti = { — Prr = ~{-)zz, Wr, = {-)zz, Vri = { — 

rj rj rj'^ rj rj rj^ 

where 

_ 1 _ w - V 

p= w = , V = -■ 

P P P 

Such reciprocal auto-transformation is the first example in the theory of integrable sys- 
tems. 

Remark 6: Twice potential two-component Harry Dym system (30) written in field 
variables (p, u) was also found in [6] (see formula (21) therein), but was not recognized. 

8. Reciprocal and Miura type transformations 

Application of the reciprocal transformation (in fact, it was given in [2], formulas 
(32-34) therein; cf. (13)) 

dyi = dti, dy_i = dt^i, d( = rjdx — rjudt_i H — ^dti (32) 

to the spectral transform (29) yields another well-known spectral transform (see for in- 
stance, [2], [27]), where "^"-dynamics is 

^CC + [A^ - - {) + = 0, (33) 

"?/" -dynamics is 

^y, = (2A + u)iJ^ - ^uci^, iJy_^ = ^(2#^ - Vci^) (34) 

and 

The compatibility conditions {4'(^i^)y-^ = ("^yjcc ii^cOy~i — ("^jz-Jcc yield the first 
positive member 

.u^ . . . 1 



yi 



25c [y + ^] ' = ^-^c [ud - -M^c] > (36) 



and the first negative member 

which also can be obtain by the limit A ^ from the generating function of commuting 
flows 

Uy = -Vc, Vy = d^[{-u - ~X)fi], 

~2 ~ ~ 1 

^CCC + '^('^^ -uX -V + -^)V(; + 2(-An^ - 0^ + -uuc)t) = 0, (37) 

where fj = ipip'^ is a squared eigenfunction and ip, ip^ are linear conjugate solutions of the 
spectral transform (33) with different asymptotics at the infinity A — > oo (see details for 
instance in [28]). 

At the same time, as it was proved in [26] (see also above section 3), the modified 
Gurevich-Zybin system is linearized by above reciprocal transformation (32). Simultane- 
ously, the tvjice 'potential two-component Harry Dym system (30) transforms into twice 
degenerate twice modified Kaup-Boussinesq system 



u 



2 



yi 



It is well-known that (36) is the Kaup-Boussinesq system (see for instance [18]). Several 
modified Kaup-Boussinesq systems were presented in [27]. The modified Kaup-Boussinesq 
system is 

_ 3 _ _ 1 _ 

^yi = ^c[2^^ + + 2^ii,d' = d(;[uiu - -m^, 

where the first Miura transformation is 

ij = ^u'^ + ul + . (39) 
The twice modified Kaup-Boussinesq system is 

^hyi = ^^([(2^1 - ui^()u2 - ^U2,(;c], U2,y^ = 9j2ni(l + u^) + U2n2,c], 
where the second Miura transformation is 

U = 2UiU2 + ^2,^ . (40) 

It was proved in [27] that the Kaup-Boussinesq system has third and fourth Miura trans- 
formations (see also [5]). Their double parametric degeneration to purely potential form 
(cf. (39), (40) with second and third equations from (35)) 

Ml = ^9^1n?7, U2 = -U(; (41) 

transforms the twice modified Kaup-Boussinesq system in form (38). 

Thus, the main result of this section is an establishment of the link of transforma- 
tions (reciprocal and Miura type) between the Gurevich-Zybin and the Kaup-Boussinesq 
hierarchies. 



9. Second Harmonic Generation 



The generation of the second harmonic wave from the red light of a ruby laser in a 
crystal of quartz in fact was the starting point of nonlinear optics. In one-dimensional case, 
for short pulses, when the group-velocity mismatch between both frequency components 
becomes important, then the process of second harmonic generation (SHG) is described 
by the complex equations (see for instance [19], [20]), which in real form are (see [20], 
formula (6) therein) 

%-i = -^c = ~2?7Ui, 2ni^j/_^ = -1 + r]u. (42) 



The corresponding spectral problem (see again [20], formula (9) therein) 



is a special reduction of the spectral transform (33), (34). In general case (37) can be 
integrated once 

1 ~2 

VVcc - + 2(A -u\-d + -)f + S{\) = 0, 

where S{\) is a polynomial function for multi-periodic solutions of the Kaup-Boussinesq 
hierarchy (see for instance [1]). Thus, the first negative member of this hierarchy has the 
constraint 

1 ^ 

where S-i 7^ is a some constant. However, S-i = in the case of the SHG system! It 
is easy to prove by direct substitution (39) and r]^ = 2rjUi from (42) (see also the first 
equation in (41)) in above equation. Thus, the SHG system (42) is the degenerate first 
negative member of the modified Kaup-Boussinesq hierarchy (see (39) and above). 
The SHG system 

{\n'^y_^)cy_^ = ^ (43) 

can be interpreted as the two-component generalization of the Sinh-Gordon equation, 
where = -u and = —77. The SHG system has three different linearizable degenera- 
tions, as well as the Sinh-Gordon equation has a parametric degeneration to the famous 
Liouville equation which is linearizable. First two degenerate limits are known (see [3] 
and [17]). These are the Liouville equation and the modified Liouville equation. The third 
such case can be obtained by (see above) differential substitutions (40) and (41). This is 
the modified Gurevich-Zybin system (31) re-calculated by the reciprocal transformation 
(32) (see formula (27) in [26] and other details in the last section 5) 

Thus, a solution of the reduced SHG system 



in an implicit form is given by (17) 

^ = ¥/ ^o(OrfC-^ I ^o(C)yo(C)rfC, 

^ = I + %(C)l/o(C)]c^C - 1 1 %(C)yo(C)^C + J 1 %(C)c^C. 

Thus, a new solution of the SHG system (43) can be found in quadratures 

Vc 

dE = —(ucc H — uc)d( — ndy-i. 
r] 

Final remark: Since the Kaup-Boussinesq system and the nonlinear Schroedinger equa- 
tion are related by invertible transformations (see for instance [1]), then their first negative 
flows are related too. Since the first negative flow to the nonlinear Schroedinger equa- 
tion is another famous Maxwell-Bloch system (in particular case the "self-induced trans- 
parency" coincides with the Maxwell-Bloch system), then the SHG system connected to 
Maxwell-Bloch system by the same transformations. Thus, the Gurevich-Zybin system is 
connected with the Maxwell-Bloch system. Similarly, a new solution of the Maxwell-Bloch 
system can be found by the same way. Since nonlinear Schroedinger equation relates to 
Heisenberg magnet by Miura type transformations, then particular case of the "Ramann 
scattering" is connecting with the Maxwell-Bloch system too. Thus, a new solution for 
the Ramann scattering can be constructed as in previous case. 



10. Open Problems 

The numerical simulation of nonlinear dynamics described by the Gurevich-Zybin 
system yields the hypothesis that a behaviour in multimode form demonstrating the 
transition from the hydrodynamic to the equilibrium kinetic state has some regular fea- 
tures (see for details [14]) possibly generated by integrable properties of corresponding 
A'"-component Gurevich-Zybin system (see (3)) 

AT 
m=l 

This problem (integrability in any sense: linearization, inverse scattering transform, bi- 
Hamiltonian formulation, etc) is open. For instance, this system written in field variables 
and t;^ {p^ = v^) has ultra- local Hamiltonian structure 

= ^, = 

where the Hamiltonian is 

m=l \m=l J 

It was proved here that the Gurevich-Zybin system has infinitely many Hamiltonian struc- 
tures. Existence of the second Hamiltonian structure is enough for an integrability. 



Introducing moments 

N 

Ak = ^^(m™)^" (44) 

m=l 

then A^— component Gurevich-Zybin system can be written as the nonlocal chain 

dtAk + d,Ak+i + kAk_id-^Ao = 0, k = 0,l,2,..., (45) 

which looks very similar as famous integrable Benney moment chain (see [4]) 

dtAk + d,Ak+i + kAk-AAo = 0, k = 0,l, 2, ... 

The Benney moment chain has infinitely many A^-component reductions (see [13]) pa- 
rameterized by functions of a single variable (A^ is an arbitrary natural integer), where 
the simplest reduction is (44). The nonlocal chain (45) has at least one simple reduction 

Ak = pu^. 

The existence of any other such reductions could be a symptom of an integrability. The in- 
tegrability of A^— component Gurevich-Zybin system and a description of other reductions 
of the nonlocal chain (45) will be under consideration anywhere else. 

11. Conclusion 

The Gurevich-Zybin system is an example of integrable systems possessing properties 
of two different classes: C— and S— integrable. This system is linearizable (see [14]) 
and has a general solution. Thus, the Gurevich-Zybin system is from a C— integrable 
class. However, this system has an infinite set of Hamiltonian structures and commuting 
fiows. Thus, the Gurevich-Zybin system is from a S— integrable class too. Moreover, this 
system has an infinite set of local Hamiltonian structures, that is unusual in the theory 
of S"— integrable systems. Moreover, all commuting fiows of the Gurevich-Zybin system 
written in a form of a Monge-Ampere equation are the same Monge-Ampere equation 
again. The difference between them is just some derivative of function p{z), which can 
be eliminated by a point transformation v = p!"{z) (see (7)). Thus, this is a beautiful 
example from mathematical point of view having physical application (see [14]). We 
repeat again: every commuting flow can be written in the unique form (23) 

where K[z) is any apriori given function, but all commuting flows will be written via 
different functions $(fc), because the point transformation like v = p"'{z) becames nonlocal 

K($(fc),.)=/i(''+')($.). 

The exceptional case is when some derivative p^^\z) is a constant: then "half of com- 
muting flows are trivial (see [22], [11]) 

when n is even, then k ^ n/2, when n is odd, then k ^ {n — l)/2. 

At past 10 years a couple of such examples of integrable systems (mixed properties of 
C— and S— integrability) was found in [22] and [12]. However, an explanation of such 
phenomenon is not exist to this moment. One possible explanation is that such systems 
are in intersection of C— and S— integrability. Thus, they accumulate properties of 
these two different classes. Moreover, we proved that the Gurevich-Zybin system is a 



degenerate member of the two-component Harry Dym hierarchy. A degeneracy becomes 
when equations can possess a parametric freedom. When some of parameters are fixed 
(to zero, for instance), then such equations become to linearizable. The simplest example 
is the famous Liouville equation 

W 

w^t = e . 

This equation is in intersection of two different integrable hierarchies. One of them is 
another famous Bonnet equation (well known in physics as the Sinh-Gordon equation) 
first introduced in a differential geometry of surfaces of a constant curvature 

Wxt = cic"" + C2e""', 

which is a member of the potential modified KdV hierarchy (spectral transform 2x2) 

1 3 

Another one is the Tzitzeica equation, well known in an affine differential geometry 

which is a member of the potential modified Sawada-Kotera hierarchy (spectral transform 
3x3) 

Wr = Wxxxxx + ^{WxxWxxx " wlw^xx - W^wl^) + W^. 

Thus, if C2 = 0, then the Liouville equation is still a member of two different integrable 
hierarchies simultaneously. This is a good symptom that such equations should be lin- 
earizable. Since above mentioned linearizable reduction of the SHG system is determined 
by 2x2 spectral transform, but the SHG system is a some reduction of another important 
three- wave interaction problem (see for instance [20]), then we can assume such systems 
like the Gurevich-Zybin system are linearizable if they are in intersection of at least two 
different integrable hierarchies (e.g. the Liouville equation is a member of two different 
hierarchies: of the KdV equation and of the Sawada-Kotera equation). A bi-Hamiltonian 
structures presented here has origin in 2x2 spectral transform. It will be interesting to 
find another Hamiltonian structures coming from 3x3 spectral problem. 

Finally, we would like to emphasize that this paper was devoted to recognition of 
a relationship between couple of remarkable systems having applications in astrophysics, 
nonlinear optics and geometry. 
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